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ABSTRACT 


The  relationships  between  the  turbulent  intensities  about  the 
mean  flow  in  the  lowest  layers  of  the  atmosphere  at  ifhite  Sands 
Missile  Range,  New  Mexico,  and  climatological  data  are  discussed  in 
terms  of  predicting  the  diffusion  of  contaminants  into  the  atmos¬ 
phere.  If  the  turbulent  fluctuations  observed  in  the  lower  portion 
of  the  boundary  layer  are  considered  ’•nrvresentative,  then  correlations 
with  climatological  summaries  may  be  used  as  a  rudimentary  estimator 
for  estimating  the  diffusion  power  of  the  atmosphere.  A  simple  em¬ 
pirical  diffusion  model  for  general  use  is  also  giver.. 
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INTRODUCTION 


Diffusion  of  aerosols  and  particulate  matter  into  the  atmosphere 
from  an  instantaneous  or  continuous  point  source  can  be  adequately 
predicted  over  short  distances  if  restrictions  are  placed  on  the  pre¬ 
diction  model.  Since  problems  of  diffusion  are  of  an  individualistic 
nature,  i.e.,  the  conditions  for  diffusion  vary  from  locale  to  locale, 
prediction  schemes  are  normally  determined  empirically  for  a  specific 
location  and  are  based  upon  the  theoretical  models  cf  Hay  and  Pasquill 
L 1 959]  or  Sutton  [1932,  1934,  1947]. 

The  development  of  a  prediction  scheme  normally  involves  detailed 
studies  of  the  area  in  question  with  respect  to  the  meteorological 
and  climatological  characteristics  and  the  use  of  tracer  materials 
and  detection  systems  to  determine  dispersion  patterns.  The  actual 
field  testing  is  normally  easy  to  accomplish;  however,  the  extraction 
of  intelligence  from  the  resulting  data  in  terms  of  the  meteorological 
situation  are  sometimes  meaningless  and  obscure.  The  major  causes  of 
a  poor  data  fit  to  theoretical  considerations  usually  stem  from  the 
fact  that  the  prediction  models  arc  based  upon  stationary  conditions 
and  homogeneous  underlying  terrain.  Since  most  locations  where  there 
is  a  diffusion  problem  are  less  than  ideal  from  a  meteorological  stand¬ 
point,  the  empirical  solutions  for  atmospheric  dispersion  of  gases 
or  matter  are  necessarily  based  upon  gross  statistics  and  regression 
techniques  and  fail  to  account  for  the  heterogeneous  effects  that  lead 
to  prediction  errors  that  can  over-  or  underestimate  by  a  factor  of  two 
to  four. 

The  scatter  about  a  regression  curve  can  be  considered  a  measure 
of  nonunifom  effects  upon  the  dispersion  of  contaminants  into  the 
atmosphere.  The  processes  determining  the  heterogeneity  of  flow  across 
a  surface  are  not  completely  understood;  consequently,  it  is  difficult 
to  compensate  for  this  in  terms  of  diffusion  theory  and  its  application. 
Thus,  timely  and  realistic  predictions  of  dispersion  corridors  and 
concentrations  of  contaminants  at  any  point  downstream  of  the  source 
must  necessarily  be  in  the  three-sigma  range  to  allow  for  the  effects 
of  nonhomogeneous  conditions. 

A  diffusion  study  for  any  locale,  then,  must  include  a  micromcteor- 
ological  study  that  considers  all  factors  of  the  mean  flow  in  steady 
and  nonsteady  state  conditions  and  over  a  wide  stability  range.  The 
vertical  fluxes  of  heat,  momentum,  and  matter  into  the  atmosphere, 
and  of  course  all  the  other  turbulent  parameters  must  be  determined  for 
the  particular  surface  over  which  the  diffusion  will  occur.  Surface 
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discontinuities  that  cause  horizontal  gradients  and  nonequi librium 
profiles  and  fluxes  great  distances  downstream  must  he  compensated 
for  in  the  analysis  of  test  data. 

PURPOSE  OF  DIFFUSION  STUDY 

The  Atmospheric  Sciences  Laboratory,  White  Sands  Missile  RangefWSMR), 
New  Mexico,  is  currently  engaged  in  a  study  of  the  meteorological 
factors  that  affect  the  diffusion  of  toxic  liquid  propellants  into  the 
atmosphere.  The  purpose  of  the  study  is  to  provide  the  National  Aero¬ 
nautics  and  Space  Administration  White  Sands  Test  Facility  (WSTF) 
with  timely  and  representative  information  concerning  possible  hazard 
corridors  and  concentration  of  toxic  vapors  in  the  vicinity  of  the 
Apollo  site. 

The  purpose  of  this  paper  is  to  present  some  preliminary  data 
concerning  the  characteristics  of  the  planetary  boundary  layer  that 
may  be  applied  to  diffusion  into  the  atmosphere.  Also  included  is 
an  empirical  diffusion  model  that  may  be  used  with  some  success  over 
a  rough,  nonhomogeneous  terrain. 

THE  SITE 

White  Sands  Missile  Range  is  located  in  South-Central  New  Mexico, 
almost  wholly  in  the  Tularosa  Basin.  The  western  boundaries  of  the  range 
extend  across  the  San  Andres  Range  of  the  southern  Rockies  into  the 
eastern  portion  of  the  Rio  Grande  Valley  and  the  Jornada  Del  Muerto. 

The  A]  olio  site  is  located  on  the  western  slopes  of  the  San  Andres 
and  is  mainly  on  the  talus  of  the  mountain  range  which  is  characterized 
by  numerous  arroyos,  rocky  soil,  and  sparse  vegetation.  The  locale 
is  heterogeneous  from  the  micrometeorological  standpoint.  Man-made 
obstructions  such  as  buildings  and  other  modifications  to  the  terrain 
contribute  to  the  overall  heterogeneity. 

TURBULENT  CHARACTERISTICS  OF  WSMR 

The  diffusive  power  of  the  atmosphere  is  directly  related  to 
the  turbulent  characteristics  of  the  planetary  boundary  layer,  and 
particularly  to  the  vertical  temperature  gradient,  the  mean  wind  speed 
and  the  direction.  The  logical  starting  point,  then,  of  any  diffusion 
study  is  to  review  the  turbulent  characteristics  of  the  atmosphere 
for  the  locale  in  question;  keeping  in  mind  of  course,  the  differences 
that  may  exist  owing  to  possible  separation  of  the  data  collection 
site  and  the  area  of  interest. 


The  investigation  of  the  turbulent  characteristics  in  the  vicinity 
of  WSMR  has  been  quite  productive  in  that  large  data  samples  have  been 
collected  and  thoroughly  analysed.  A  goodly  portion  of  the  data  are 
applicable  to  the  Apollo  side  of  the  mountains,  even  though  they  were 
collected  in  the  Tularosa  basin  portion  of  WSMR. 

It  has  been  shown  by  Monin-Obukhov  [1954],  Panofsky  [1961,  1963] 
Ellison  [1957]  and  others  that  atmospheric  turbulence  obeys  the  laws 
of  dynamic  similarity  of  flows.  In  the  surface  boundary  layer  and 
to  lesser  extent  in  the  planetary  boundary  layer  proper,  the  wind  and 
temperature  profiles  can  be  shown  to  be  given  by,  respectively, 

\ 

r  • s  Ri  VSi  (i) 

and 

rRRi  W2  (2) 

where  z/L  is  the  Monin-Obukhov  scaling  ratio,  Ri  is  the  gradient 
Richardson  number,  S  is  a  dimensionless  wind  shear,  R  is  a  dimension¬ 
less  lapse  rate  and  Kh/Km  is  the  ratio  of  the  exchange  coefficients 
for  heat  and  momentum,  batchelor  [1953]  has  shown  that  the  Richardson 
number  is  the  basic  parameter  of  similarity  for  the  boundary  layer,  and 
it  has  been  demonstrated  by  Hansen  [1966],  and  Swanson  and  Cramer 
[1965]  that  the  turbulent  characteristics  of  the  boundary  layer  are 
directly  related  to  this  basic  stability  parameter.  The  Richardson 
number  is  usually  expressed  as 


Ri 


JL  3Q/3z 
F  (aV/'O 


(3) 


where  g  is  the  acceleration  of  gravity,  7  is  potential  temperature, 
V  is  the  mean  horizontal  wind  speed  and  z  is  height.  From  *he  simi¬ 
larity  theory  it  may  be  shown  that  the  Richardson  number  is  also  a 
function  of  height  and  its  vertical  distribution  governed  by 
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where  (Ri)  denotes  the  height  derivative  of  Hi.  An  investigation 
of  the  distribution  of  Ri  with  time  as  well  as  height  by  Hansen  and 
Lang  [1966]  reveals  that  eq.  (4)  adequately  describes  the  height 
distribution  of  Ri  to  the  top  of  the  surface  boundary  layer  which, 
depending  upon  stability,  is  some  10  to  40  netors  above  the  surface. 
This  is  shown  in  Figures  1  and  2  for  both  homogeneous  and  nonhono- 
geneous  flow  conditions.  Figures  1  and  2  reveal  that  critical  values 
of  Ri  exist,  especially  in  the  unstable  regime.  It  has  been  found 
by  Dyer  [1965]  and  Hansen  [1967]  that  critical  Ri  is  in  the  range 
-0.70  >  Ri  >  -1.2.  In  the  stable  regime,  critical  Iti  is  thought  to 
lie  between  Ri  «  0.20  and  Ri  =  0.40.  Outside  these  limits  arc  regimes 
of  local  r.onstationarity  known  as  the  "undulant"  repine*  for  very 
stable  conditions  and  the  "windless  convection"  regime  for  the  ex¬ 
tremely  unstable  case. 

It  was  also  found  that  temporally  the  distribution  of  Ri  with 
height  was  as  depicted  in  Figure  3.  It  is  seen  that  the  Richardson 
number  varies  with  the  energy  balance  of  the  boundary  layer,  which 
is  best  shown  in  Figure  4.  From  Figures  3  and  4  it  is  seen  that  the 
atmosphere  is  almost  wholly  thermally  stratified  at  all  times  and  that 
neutral  conditions  exist  only  for  very  short  periods  near  sunrise  and 
sunset.  It  is  interesting  to  note  that  the  Richardson  number  at 
3  and  6  meters  is  almost  perfectly  correlated  with  the  incoming  short 
wave  and  net  radiation  curves,  indicating  that  the  stability  of  the 
atmosphere  at  these  heights  reflectr'steady-state  conditions  much 
better  than  some  other  reference  height,  say  1.5  meters. 

If  the  basic  similarity  parameter  is  taken  to  he  the  Richardson 
number,  and  if  it  is  taken  to  be  an  independent  variable  in  the 
boundary  layer,  cqs.  (1)  and  (2)  may  be  integrated  to  find  the  wind 
and  temperatuie  profiles.  Substituting  in  and  integrating  cq.  (1) 
yields 


V  “  T~  lln  f“  (I/L)1  (5) 

o 

which  is  known  as  the  log-linear  wind  prof ilc,  where  u#  is  a  stress 
term,  k  is  Karnan's  constant,  iQ  is  the  roughness  length  and  ii»(t/l.) 
is  a  universal  function.  A  similar  expression  may  be  written  for 
the  temperature  profile,  although  not  in  such  a  universal  form. 

The  turbulent  characteristics  of  the  boundary  layer  nay  then  be 
expressed  as  deviations  about  the  mean  profile  as  functions  of  Pi. 
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FIG  3  THE  RICHARDSON  DUMBER  AS  A  FUNCTION  OF  HEIGHT  AND  TIME 
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FIS  4  THE  RICHARDSON  NR.  AT  SELECTED  HEIGHTS  VERSUS  NET  AND  INC0MSN6  SHORT 
WAVE  RADIATION 


Swanson  an<J  Cramer  [1965],  Hansen  [1966]  and  Hansen  and  Lang  4 

[1967]  have  shown  that  the  turbulent  wind  and  temperature  fluctuations 
about  the  means  are  directly  related  to  the  stability  regime.  It  was  > 

found  that  the  orthogonal  components  of  the  horizontal  wind  were  i 

functions  of  stability  and  to  a  lesser  extent  functions  of  surface  | 

roughness.  The  turbulent  intensity  was  also  found  to  he  height  dc-  \ 

pendent  in  the  form 


Gx,  Gy  -  Czp  (6) 

where  Gx  and  Gy  are  defined  as 

Gx  ■  ou  /V  (7) 

o 

Gy  •  ov  /V  £  oAjj  (8) 


where  ou*  and  ov*  are  the  standard  deviations  of  the  longitudinal 
and  lateral  eddy  velocities  and  oAr  is  the  standard  deviation  about 
the  mean  wind  direction  in  radians.  The  constant  C  is  a  function 
of  surface  roughness  while  the  exponent  p  is  related  to  the  exponent 
of  the  conjugate  power  law  wind  profile  and  is  a  function  of  Ri. 

Hansen  and  Lang's  [1967]  results  are  shown  in  Figures  5  through 
8.  It  will  he  noted  that  an  inflection  occurs  in  the  oA  versus  z 
regression  of  Figure  8  for  the  stable  case.  This  is  attributed  to 
a  surface  roughness  discontinuity  upstream  of  the  research  tower  as 
reported  by  Hansen  and  Hansen  [1965]  .  The  inflection  in  the  curve 
is  the  result  of  excessive  surface  stress  downstream  of  the  dis¬ 
continuity  and  the  lack  of  vertical  transport  of  momentum  in  stable 
flow  conditions  through  the  internal  momentum  boundary. 

Turbulent  temperature  fluctuations  over  a  nonhonogeneous  surface 
are  shown  in  Figure  9  as  a  function  of  the  Richardson  number.  It 
is  seen  that  in  the  region  below  the  internal  momentum  boundary  the 
quantity  oT/30/3z  docs  not  exhibit  the  predicted  characteristics 
as  indicated  by  Tsvang  [1960]  and  Lumley  and  Panofsky  [1964].  Above 
the  interface  height  of  8  to  10  meters,  the  scaling  temperature 
T*(T*%  aT/30/3z)  has  the  predictable  characteristic  of  being  nearly 
unity  in  unstable  flow  and  independent  of  height  in  stable  conditions. 
The  dashed  portions  of  this  curve  indicate  local  nonstationav ity 
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FIGURE  7 


HEIGHT  DEPENDENCE  OF  THE  LATERAL  INTENSITY  OF 
TURBULENCE  DURING  UNSTABLE  STRATIFICATION 
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effects  upon  the  experimental  data.  On  the'  stable  side,  it  appears 

J  that  undulance  has  the  same  effects  as  the  nonequilhrium  'low  beneath 

the  boundary. 

t  TERRAIN  EFFECTS  ON  THE  MEAN  WIND 

r 

j  The  terrain  of  WSMR  and  particularly  the  proximity  os'  a  mountain 

range  to  the  experimental  site  has  a  profound  effect  upon  the  annual 
wind  regime  at  WSMR.  Hansen  and  Lang  [1966b,  1966c]  four,  that  the 
annual  wind  regime  could  be  expressed  in  terms  of  three  s  asonal 
regimes.  These  consisted  of  a  winter  regime,  a  two-part  .pring 
regime  and  a  combined  summer-fall  regime.  The  winter  ami  spring 
seasons  are  characterized  by  generally  strong  westerly  fl  w  with 
the  summer-fall  regime  being  generally  southeasterly.  All  seasons 
reflect  a  significant  percentage  of  drainage  winds  from  t'.c  mountain 
slopes  during  the  nocturnal  hours.  The  percent  frequency  of  occurrence 
of  wind  direction  for  all  seasons  is  shown  in  Figures  10  through 
13. 


Other  terrain  effects  such  as  the  roughness  length  have  been 
investigated.  Hansen  [1967]  has  found  that  the  roughness  length  of 
the  terrain  about  the  research  tower  varies  with  wind  direction  and 
with  the  season.  The  seasonal  change  in  roughness  is  shown  in  Figure 
14.  It  will  be  noted  that  the  curve  approximates  the  growing  season 
in  that  the  largest  values  of  zQ  are  found  in  July  and  August  when 
the  vegetive  canopy  is  most  dense. 

AN  EMPIRICAL  DIFFUSION  PREDICTION  MODEL 

Experimental  evidence  shows  that  three  meteorological  parameters 
should  be  considered  in  characterizing  the  rate  of  small  scale  at¬ 
mospheric  diffusion: 

(1)  The  mean  wind  speed  V  as  an  indicator  of  the  down¬ 
stream  travel  of  the  plume  emitted  from  a  continuous  source. 

(2)  The  standard  deviation  of  wind  direction  o A  as  an 
indicator  of  the  lateral  rate  of  mixing. 

(3)  The  temperature  gradient  &T  as  an  indicator  of  the 
vertical  rate  of  mixing. 

Other  quantities  of  interest  are  the  peak  or  maximum  concentration 
Cp  at  some  distance  x  downstream  of  the  source  and  the  related  source 
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FIGURE  II  PERCENT  FREQUENCY  OF  OCCURRENCE  OF  WIND  DIRECTION,  SPRING  (PART  I 


FIGURE  12  PERCENT  FREQUENCY  OF  OCCURRENCE  OF  WIND  DIRECTION,  SPRING  (PART  2) 
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FIGURE  !3  PERCENT  FREQUENCY  OF  OCCURRENCE  OF  WIND  DIRECTION,  SUMMER  -  FALL 


MEXMCO 


strength  Q.  In  the  development  of  a  regression  equation  for  pre¬ 
dicting  diffusion,  a  normalized  peak  concentration  Cp/Q  at  distance 
x  is  generally  the  goal. 

The  regression  equation  usually  takes  the  form 

Cp/Q  »  Kxa  V  b  oAC(AT*t)d  (0) 

where  K,  a,  b,  c,  d  are  parameters  of  fit  and  t  is  a  constant  added 
to  AT  tc  avoid  raising  a  negative  number  to  a  power. 

Application  of  linear  multiple  regression  techniques  requires 
taking  the  logarithm  of  eq  (9)  and  applying  standard  least-squares 
analysis  to  the  following: 


log  C^Q  ■  log  K  ♦  a  log  x  ♦  b  log  V  ♦  c  log  oA  ♦  d  log(AT+t)  (10) 

since  the  normalized  peak  concentration  could  be  linear  with  respect 
to  the  meteorological  variables  but  not  linear  in  distances. 

A  general  purpose  form  of  equation  (9)  for  use  in  nonhomogcncous 
terrain  has  been  developed  by  Miller  and  Miller  [1963]  in  the  form 


Cp/Q  *  0.000175  x"1,95(AT*10)4,95  (11) 

where  Cp  is  the  peak  concentration  in  gms  n”0  at  a  height  of  about 
1.5  meters  above  the  surface  at  a  downwind  travel  distance,  x,  in 
meters,  Q  is  the  source  strength  in  gn  see”*  and  AT  is  the  temperature 
difference  in  deg.  F  between  16.5  and  1.8  meters  above  the  ground. 

Eq.  (11)  was  developed  from  experimental  data  simulating  ground-level 
continuous  point  sources. 

It  will  be  noted  that  the  mean  wind  speed  V  and  the  standard 
deviation  of  wind  direction  are  not  included  in  eq.  (11).  It  was 
found  by  Haugen  and  Taylor  [1963]  and  Miller  and  Miller  [1963]  that 
while  7  is  correlated  with  Cp/Q,  its  contribution  to  the  predicta¬ 
bility  after  the  temperature  differential  is  included  is  not  signifi¬ 
cant.  The  oA  term  was  dropped  owing  to  the  fact  that  i  satisfactory 


method  of  obtaining  oA  from  an  analogue  recorder  was  available. 
Hstimates  for  oA  will  be  introduced  below. 

For  direct  application,  cq  (11)  must  be  converted  to  units  which 
may  be  used  conveniently.  The  conversion  is  relatively  simple  for 
contaminants  emitted  as  a  ground-level  continuous  point  source.  In 
the  case  of  a  gaseous  contaminant,  Cp/Q  should  be  known  in  terns 
of  parts  per  million  of  the  contaminant  per  pound  per  minute  released 
into  the  atmosphere,  a  value  that  is  dependent  upon  the  diffusive 
power  of  the  atmosphere  and  the  physical  characteristics  of  the  con¬ 
taminant.  Under  standard  conditions,  1  ppm  lb"1  min_equals  the  gran 
molecular  weight  of  the  gas  times  5. 89  x  10"^  sec  m"^._  To  obtain 
parts  per  million  per  pound  per  minute,  Cp/o  in  sec  m“  ,  is  multi¬ 
plied  by  1.70  x  1CP  and  divided  by  the  gram  molecular  weight. 

In  terms  of  travel  distance  for  normalized  peak  concentrations 
eq  (11)  nay  be  restated  as 


x(in  feet)  =  2.62(Cp/q)'°'513  [AT-10]2*55  (12) 

and 

x(in  feet)  «=  2.20  (Cp/Q)-0*513  [AT«-10]2*53  (13) 

Lqs.  (12)  and  (13)  arc  valid  for  two  well-known  propellants,  nitrogen 
tetroxidc  and  unsynmetr ical  dimethylhydrazinc.  For  other  fuels  or 
oxidizers,  the  constants  of  eqs.  (12)  and  (13)  must  be  modified  accord¬ 
ing  to  the  scheme  outlined  above. 

Travel  tine  of  the  contaminant  to  a  particular  location  can 
be  determined  from  the  mean  wind  speed.  The  width  of  the  hazard 
corridor  can  be  estimated  from  a  wind  direction  trace  by 

W  «  C  R  (14) 

where  W  is  the  width  of  the  corridor,  !l  is  the  range  of  the  wind 
direction,  and  C  is  an  empirical  constant  determined  for  a  locale  and 
a  particular  recorder  speed,  l.q,  (14)  yields  equivalent  results  obtained 
by  Haugen  and  Taylor  [1963]  who  showed  that  the  width  of  the  hazard 
corridor  could  be  estimated  from  K  »  8  oA  with  95*  confidence  that 
this  value  would  not  be  exceeded. 


DISCUSSION  AND  SUMMARY 

•j 

If  the  turbulent  characteristics  of  the  boundary  layer  as  observed 
at  KSMR  can  be  considered  typical  of  an  area  that  includes  the  WSTF 
(Apollo)  locale,  some  insight  into  the  diffusion  of  a  contaminant 
into  the  atmosphere  is  provided.  If  the  KSMR  turbulence  intensities 
are  typical,  then  climatological  type  data  at  WSTF  may  be  used  to 
investigate  the  diffusive  power  of  the  atmosphere,  alleviating  any 
need  for  elaborate  multilevel  tower  systemswhich  are  a  necessity 
at  KSMR. 

The  only  exception  to  this  approach  would  be  the  drainage  wind 
pattern  during  the  nocturnal  hours.  It  is  clear  that  the  drainage 
pattern  at  WSTF  would  not  be  the  same  as  it  is  at  WSMR.  This  is 
because  WSSfR  and  WSTF  are  on  opposite  slopes  of  the  mountain  range. 

It  appears  likely  that  the  nocturnal  drainage  on  the  western  slopes 
of  the  San  Andres  would  be  almost  a  mirror  image  of  the  eastern  slope 
drainage  pattern  with  few  exceptions. 

The  use  of  the  Miller  and  Miller  [1963]  diffusion  equation 
for  predicting  concentrations  and  hazard  corridors  for  WSTF  is  perhaps 
not  completely  justified  but  is  a  useful  tool  for  estimating  the 
dispersion  of  contaminants  in  the  atmosphere.  The  development  of  a 
more  precise  empirical  model  for  WSTF  would  entail  a  very  complex 
and  costly  experiment. 

This  analysis  indicates  that  the  turbulent  characteristics  of  the 
atmosphere  as  observed  at  White  Sands  arc  representative  of  larger 
areas.  Under  these  conditions  then,  the  diffusive  power  of  the 
atmosphere  at  the  Apollo  site  may  be  estimated  from  rudimentary  cli¬ 
matological  data  used  in  conjunction  with  the  KSMR  data. 
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The  relationships  between  the  turbulent  intensities  about  the  mean  flow 
in  the  lowest  layers  of  the  atmosphere  at  White  Sands  Missile  Range,  New  Mexico, 
and  climatological  data  are  discussed  in  terms  of  predicting  the  diffusion  of 
contaminants  into  the  atmosphere.  If  the  turbulent  fluctuations  observed  in  the 
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tions  with  climatological  summaries  may  be  used  as  a  rudimentary  estimator  for 
estimating  the  diffusion  power  of  the  atmosphere.  A  simple  empirical  diffusion 
model  for  general  use  is  also  given.  , 
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